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ABSTRACT 

We present an analysis of the color-magnitude- velocity dispersion relation for a sample of 39320 early- 
type galaxies within the Sloan Digital Sky Survey. We demonstrate that the color-magnitude relation is 
entirely a consequence of the fact that both the luminosities and colors of these galaxies are correlated 
with stellar velocity dispersions. Previous studies of the color-magnitude relation over a range of redshifts 
suggest that the luminosity of an early-type galaxy is an indicator of its metallicity, whereas residuals 
in color from the relation are indicators of the luminosity-weighted age of its stars. We show that this, 
when combined with our finding that velocity dispersion plays a crucial role, has a number of interesting 
implications. First, galaxies with large velocity dispersions tend to be older (i.e., they scatter redward of 
the color-magnitude relation). Similarly, galaxies with large dynamical mass estimates also tend to be 
older. In addition, at fixed luminosity, galaxies which are smaller, or have larger velocity dispersions, or 
are more massive, tend to be older. Second, models in which galaxies with the largest velocity dispersions 
are also the most metal poor are difficult to reconcile with our data. However, at fixed velocity dispersion, 
galaxies have a range of ages and metallicities: the older galaxies have smaller metallicities, and vice- 
versa. Finally, a plot of velocity dispersion versus luminosity can be used as an age indicator: lines of 
constant age run parallel to the correlation between velocity dispersion and luminosity. 

Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: photometry — galaxies: stellar content 



1. INTRODUCTION 

It has long been known that early-type galaxies tend to 
be the reddest galaxies, and that the more luminous the 
galaxy, the redder its color (e.g. Sandage & Viswanathan 
1978a,b; Bower, Lucey & EUis 1992a,b). The tightness of 
the correlation between color and magnitude has been used 
to constrain models of how early-type galaxies formed. 
However, a luminous galaxy may appear red either be- 
cause its stars are older, or because, although its stars are 
younger, they are more metal rich. This has complicated 
the constraints one can place on galaxy formation models: 
should the models produce luminous metal-rich galaxies, 
or luminous old galaxies? To illustrate. Appendix A shows 
the color-magnitude relation associated with the recent 
stellar population synthesis models of Bruzual & Chariot 
(2003). 

The colors and luminosities of older galaxies are ex- 
pected to evolve more slowly than those of younger ones. 
Therefore, if age is changing along the color-magnitude se- 
quence (e.g. if the more luminous galaxies are older), one 
would expect the differential evolution of the older and 
younger populations along the sequence to manifest as a 
change in the slope of the color-magnitude relation with 
redshift. Recent measurements have shown that the rela- 
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tion was already in place at redshifts of order unity, and 
that its slope appears to be unchanged from its value lo- 
cally (e.g. Kodama et al. 1998; Blakeslee et al. 2003). 
If the high redshift population does indeed represent the 
local low redshift population in its youth, then the fact 
that the slope has evolved little argues against a large 
age-spread along the relation. In this case, the color- 
magnitude relation is caused primarily by a correlation 
between metallicity and luminosity (Kodama et al. 1998). 
(This argument becomes weaker if the redshift at which 
the population formed the bulk of its stars is large.) 

In this framework, the scatter around the color-magnitude 
relation is usually attributed to the effects of age. It is 
not entirely obvious, however, that this can be correct. 
This is because differential evolution must make the scat- 
ter around the mean relation larger at higher lookback 
times. Therefore, the fact that the color-magnitude re- 
lation is well defined at redshifts of order unity can be 
translated into a constraint on the mix of ages present at 
redshift zero. For this reason, it would be interesting to 
quantify how the scatter around the mean relation evolves. 
(Blakeslee et al. 2003 find little evolution in the slope and 
scatter of the relation, using cluster early-types out to red- 
shifts of order unity.) 

It is well known that luminosity also correlates with ve- 
locity dispersion (Poveda 1961; Faber & Jackson 1976), 
so it is natural to ask if velocity dispersion is also tightly 
coupled to metallicity. The tightness of the correlation 
between Mg2 and a (e.g. Bernardi et al. 1998; CoUess et 
al. 1999), is thought to be a consequence of variations in 
age and metallicity which conspire to keep the observed 
correlation tight (e.g. Kuntschner et al. 2001). In this 
view, old galaxies with a given a tend to be metal poor, 
whereas younger galaxies with the same a are metal rich 
(e.g. Worthey 1994; Trager et al. 2000b). 
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Thus, recent work suggests that residuals from the color- 
magnitude relation are age indicators, whereas studies of 
chemical abundances suggest that velocity dispersion is a 
combination of both age and metallicity. One of the goals 
of the present work is to see if these conclusions are con- 
sistent with one-another. We do this by studying the joint 
distribution of color, luminosity and velocity dispersion 
in a sample of 39320 early-type galaxies drawn from the 
Sloan Digital Sky Survey (York et al. 2000; Stoughton 
et al. 2002; Abazajian et al. 2003) database (hereafter 
SDSS). 

Section 2 describes our sample. Section 3 presents the 
color-magnitude and color-cr relations, and demonstrates 
that color-fT and magnitude-cr are the primary correlations. 
A simple model is introduced which illustrates clearly what 
our measurements imply for the relations between age, 
metallicity and velocity dispersion. The mathematics as- 
sociated with this model is in Appendix B. Section 4 uses 
the model to derive a number of consequences of the ob- 
served color-magnitude-(T relations, and summarizes our 
findings. 

Throughout, we assume that Hq — 70 km s^^ Mpc^^ in 
a universe with Qq = 0.3 which is spatially flat. 

2. THE SDSS EARLY- TYPE GALAXY SAMPLE 

For our analysis we used galaxies selected from the 
Sloan Digital Sky Survey (SDSS) database. See York et 
al. (2000) for a technical summary of the SDSS project; 
Stoughton et al. (2002) for a description of the Early Data 
Release; Abazajian et al. (2003) et al. for a description 
of DRl, the First Data Release; Gunn et al. (1998) for 
details about the camera; Fukugita et al. (1996), Hogg et 
al. (2001) and Smith et al. (2002) for details of the photo- 
metric system and calibration; Lupton et al. (2001) for a 
discussion of the photometric data reduction pipeline; Pier 
et al. (2003) for the astrometric calibrations; Blanton et 
al. (2003) for details of the tiling algorithm; Strauss et al. 
(2002) and Eisenstein et al. (2001) for details of the target 
selection. 

Our early-type galaxy sample is selected similarly to 
that of Bernardi et al. (2003a), with some minor changes: 
We selected all objects targeted as galaxies and having 
de-reddened Petrosian apparent magnitude 14.5 < rpct < 
17.75. To extract a sample of early-type galaxies we then 
chose the subset with the spectroscopic parameter eel ass 
< (eclass classifies the spectral type based on a Princi- 
pal Component Analysis), and the photometric parameter 
fracDeVj. > 0.8. (The parameter fracDev is a seeing- 
corrected indicator of morphology. It is obtained by tak- 
ing the best fit exponential and de Vaucouleurs fits to the 
surface brightness profile, finding the linear combination 
of the two that best-fits the image, and storing the fraction 
contributed by the de Vaucouleurs fit.) We removed galax- 
ies with problems in the spectra (using the zStatus and 
zWarning flags). From this subsample, we finally chose 
those objects for which the spectroscopic pipeline had 
measured velocity dispersions (meaning that the signal- 
to-noise ratio in pixels between the restframe wavelengths 
4200A and 5800A is S/N > 10). These selection criteria 
produced a sample of 39320 objects with photometric pa- 
rameters output by version V5.4 of the SDSS photometric 
pipeline and V.23 reductions of the spectroscopic pipeline 
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Fig. 1. — Luminosity correlates with both velocity dis- 
persion (bottom) and size (middle), but the tightest cor- 
relation is with Rra^ oc mass (top). Dashed lines show 
the correlations summarized in Table 1 . Luminosities have 
been corrected for evolution by 0.85z (see text for details). 



(i.e. v4_10_2 of idlspec2d and v5_9_4 of spectrold). 

The main quantities we use in this paper are magni- 
tudes, sizes, colors, redshifts and velocity dispersions; the 
first three are output from the SDSS photometric pipeline, 
and the final two from the SDSS spectroscopic pipeline. 
The photometry is accurate to 0.02 mags, and the veloc- 
ity dispersions to about 8% (see Bernardi et al. 2003a). 
We work with de Vaucouleurs (1948) magnitudes and sizes, 
and SDSS model colors throughout: these are colors mea- 
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Table 1 

Maximum-likelihood estimates, in the SDSS r band, of the joint distribution of luminosities, sizes and velocity dispersions. 
The mean values of the variables at redshift z, M* — Qz, i?*, 14, and the coefficients of the various pairwise correlations 
between the variables shown. TIk^ coc^fficic^nt which differs most from Table 1 of Bc^rnardi (^t al. (2003b) is i?*. 



Band M* 




R* 






aw 


Crm 




^RV 


Q 


r -21.025 


0.841 


0.354 


0.263 


2.220 


0.111 


-0.900 


-0.774 


0.550 


0.85 



sured within an aperture which scales with the de Vau- 
couleurs half-light radius in the r-band. 

The most important reason for analyzing a new sam- 
ple (other than size) is that the old photometric reduc- 
tions output by the SDSS pipeline photo were incorrect 
(sec documentation on DR2, the Second Data Release). 
A comparison of the properties of objects for which old 
and new photometric reductions are available shows that 
the new corrected photometry has made most magnitudes 
slightly fainter (^ 0.13 mags), and most half-light radii 
smaller (~10%). It is customary to report velocity disper- 
sions at some fraction (we use 1 /8) of the half-light radius. 
The half-light radii of the galaxies which enter our sample 
are typically about 2 arcsecs, approximately independent 
of redshift, whereas the SDSS fiber used to measure the 
spectrum from which the velocity dispersion is estimated 
has a diameter of 3 arcsec. Since the velocity dispersions of 
early-type galaxies are known to increase towards the cen- 
ter (following J0rgenscn et al. 1995 we assume the scaling 
is oc (r/re)°'°^, where is the half-light radius), all mea- 
sured velocity dispersions are 'aperture corrected', and it 
is these which are usually reported (the mean correction is 
7%). Although the measured velocity dispersions have not 
changed, the new photometry has changed the half-light 
radius. Hence, aperture corrected velocity dispersions dif- 
fer from those associated with the old photometry (i.e., 
those reported in Bernardi et al. 2003a) by less than one 
percent (the new values are larger by a factor of l.l" ""'). 

Bernardi et al. (2003b) showed that the luminosities 
in their sample (~ 9000 galaxies with z < 0.3) evolve: 
M^{z) = M«(0) — 0.85^;, where M^,{z) denotes the mean 
absolute magnitude at redshift z. We find similar evolu- 
tion in the new sample: the main difference is that the 
new photometric reductions make fainter by about 
0.125 mags. The evolution is consistent with that of a 
passively aging population. Thus, in the analysis which 
follows, we will be careful to separate trends with red- 
shift that are due to the magnitude limit of the sample, 
from trends that are due to evolution. Many of the fol- 



Table 2 

Maximum-likelihood estimates of the joint distribution of 
color, r-band magnitude and velocity dispersion and its 
evolution. At redshift z, the mean values are C* — Pz, and 

the covariances are ((C — C*)(M — M*)) = acc^^MM S,cm, 
and similarly for i^cv {^* and 14 are given in Table 1). 
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g-r 


0.736 


0.057 


-0.361 


0.516 


0.30 



lowing figures present measurements made in a scries of 
narrow redshift bins: 0.02 < z < 0.07, 0.07 < z < 0.09, 
0.09 < z < 0.12, 0.12 < z < 0.15, and 0.15 < z < 0.2. In 
all of these figures, symbols show the median value of the 
variable on the y-axis as a function of the observable on 
the x-axis, error bars show the rms error of the median, 
and hashed regions indicate where the central 68% of the 
sample in each redshift bin lies. 

Figure 1 displays correlations between luminosity, size 
Rr, and velocity dispersion a in this sample. (The half- 
light radius in kpc, Rr, is defined similarly to the quantity 
called Ro by Bernardi et al. 2003a, for the r-band.) The 
dashed lines show fits to these correlations, summarized 
in Table 1. In all cases, luminosities have been corrected 
for evolution by 0.85^. In the Table, M» and ctm are the 
mean and rms values of a Gaussian fit to the evolution- 
corrected r-band luminosity function, Q quantifies the rate 
of luminosity evolution, and iZ* and K are the values of 
log]^Q(size//i~"'^kpc) and logj^Q (velocity dispersion/km s"-'^) 
for an galaxy. In general, these values are all qualita- 
tively similar to those reported in Bernardi et al. (2003b), 
but there are quantitative differences, most noticably in 
the mean size R^,, which is smaller by about ten percent. 
Table 1 also reports values of a number of pairwise correla- 
tions in the following format. Given a pair of observables 
X andy, {Y-Y,\X)/aYY = ^xy{X - X,)/axx. There- 
suits of a similar analysis of the correlations between color, 
magnitude and velocity dispersion is given in Table 2. 

3. COLOR-MAGNITUDE AND COLOR-a RELATIONS 

Figure 2 shows the color-magnitude relation in our sam- 
ple. The top series of panels shows the raw measurement, 
and bottom panels show the result of accounting for evolu- 
tion by adding 0.852; mags to the magnitudes, and adding 
0.3z mags to the colors (see Table 2). On both tiers, the 
left-most panel shows the color-magnitude relation in the 
full sample, and the other two panels show the relation in 
small bins in velocity dispersion. The figure makes two im- 
portant points: (i) at fixed velocity dispersion, there is no 
relation between color and magnitude, and (ii) the slope 
of the relation in the leftmost panel is approximately the 
same in all redshift bins. 

We also find no evidence that the scatter around the 
mean relations changes with redshift, but placing a pre- 
cise limit of the evolution is difficult, because the expected 
signal is small, so the answer is sensitive to uncertainties 
in the SDSS photometry (0.02 mags) and the fc-correction. 
Specifically, because we fc-correct following the procedure 
described in Bernardi et al. (2003a), all galaxies at the 
same redshift are assigned the same fc-corrcction. This 
leads to a small amount of scatter (less than 0.02 mags) 
in our absolute magnitudes and colors which may well de- 
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Fig. 2. — Color-magnitude relation in the sample. Top series of panels show the raw measurement, and bottom panels 
show the result of making magnitudes fainter by 0.85z mag and colors redder by 0.3z mag to account for evolution. Both 
on the top and bottom, the left-most panel shows the color-magnitude relation in the full sample, and the other two panels 
show the relation in small bins in velocity dispersion. Clearly, at fixed velocity dispersion, there is no relation between 
color and magnitude. 



pend on redshift. 

Figure 3 shows the color-cr relation in the sample. The 
format is the same as for the previous figure: top panels 
show the raw measurement, and bottom panels show the 
result of accounting for evolution by adding 0.32: mags to 
the colors, as was required to model the evolution of the 
color magnitude relation. The left-most panels show the 
color-<T relation in the full sample, and the other two pan- 
els show the relation for galaxies which are restricted to 
a narrow range in magnitude. Clearly, at fixed redshift, 
the correlation between color and cr is the same for all 
magnitude bins. 

Comparison of the color-cr relation in the different red- 
shift bins suggests that the relation is steepening slightly 
with redshift: galaxies with small velocity dispersions are 
evolving slightly more rapidly. Because we are no longer 
making measurements at fixed magnitude, one may won- 
der if the slope of the color-cr relation is affected by the 
magnitude-limited selection. Appendix B develops a sim- 
ple model which demonstrates that this is not the case. We 
have also checked that the slope change does not depend 
strongly on our choice of fc-correction. With the current 
sample size, this steeping is barely significant: it will be in- 
teresting to see if this steepening persists when the sample 
is larger. 

Together, Figures 2 and 3 indicate that the color- 



magnitude relation is entirely a consequence of the color-cr 
and magnitude-cr relations. In particular, the joint distri- 
bution of color, magnitude and velocity dispersion is: 

p(color, mag, cr) = p(cr) p(mag|cr) p(color|cr). (1) 

Appendix B demonstrates that if this expression is correct, 
then residuals from the magnitude-cr relation should not 
correlate with residuals from the color-cr relation. This 
is shown in Figure 4: M — (M|cr) does not correlate 
with color — (colorjcr), but a — {cr\M) does correlate with 
color — (colorjM). Moreover, residuals from the color- 
cr relation should not correlate with magnitude, whereas 
residuals from the color-magnitude relation should corre- 
late with velocity dispersion. The two panels in Figure 5 
show that this is indeed the case. Note that the slope 
of the correlation between color — (colorjM) and a is not 
as steep as for the color-cr relation itself. In effect, this 
is because (color |Af) correlates with cr (because M itself 
correlates with a). 

The bottom panel of Figure 5 shows that the correlation 
between color-magnitude residuals (i.e., age) and velocity 
dispersion steepens systematically with redshift. Although 
it is tempting to conclude that this is evidence of differen- 
tial evolution, the analysis in Appendix B shows that it is, 
in fact, a selection effect (see text between equations B9 
and BIO). 
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Fig. 3. — Color-cr relation in the sample. Top series of panels show the raw measurement, and bottom panels show the 
result of accounting for evolution by adding 0.3z mags to the colors, as was required to model the evolution of the color 
magnitude relation. Both on the top and bottom, the left-most panel shows the color-cr relation in the full sample, and 
the other two panels show the relation in small bins in magnitude. Clearly, at fixed redshift, the correlation between color 
and a is the same, whatever the magnitude. 



4. INTERPRETATION AND DISCUSSION 
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Fig. 4. — Residuals from the color-magnitude relation cor- 
relate with residuals from the u-magnitude relation (left) 
whereas residuals from the color-cr relations do not corre- 
late with residuals from the magnitude-a relation (right). 
Dashed lines show the relations from our model, equa- 
tion (Bl), obtained by inserting the values from Tables 1 
and 2 in equations (B7) and (B6), respectively. 



Suppose that the color-magnitude relation is a conse- 
quence of the fact that metallicity increases with luminos- 
ity along the sequence, and that the galaxies which scatter 
redward of the relation are older. Then the correlation in 
Figure 5 suggests that galaxies with large velocity disper- 
sions are, on average, older. 

The correlation in the panel on the left of Figure 4 im- 
plies that older galaxies tend to have larger velocity dis- 
persions than their luminosities indicate. A correlation 
between age and the residuals from the cr-luminosity rela- 
tion was found by Forbes & Ponman (1999), but their age 
estimates came from comparing the spectra of these ob- 
jects with stellar population synthesis models. Our analy- 
sis indicates that age estimates from the color-magnitude 
relation lead to the same conclusion. (Note that the analy- 
sis in Appendix B suggests that age is more closely related 
to the residuals in cr with respect to the mean (cr|L) rela- 
tion, than to the residuals in magnitude with respect to the 
{L\a) relation. That is, the residuals studied by Forbes & 
Ponman (1999) should be more closely related to age than 
those studied more recently by Proctor et al. 2004.) 

Figure 1 shows that luminosity is tightly correlated with 
mass oc R^cr^ . So it is interesting to study how residuals 
from the color-magnitude relation correlate with residuals 
from the mass-luminosity relation. The top-most panel of 
Figure 6 shows that galaxies which scatter blueward from 



6 



o 



0.10 



0.05 



0.00 



-0.05 



-0.10 




-20 -21 -22 -23 
[mag] 




2.0 2.2 2.4 2.6 
Log,o a 



Fig. 5. — Residuals from the color-cr relations do not 
correlate with magnitude (upper panel) , whereas residuals 
from the color-magnitude relation do correlate with ve- 
locity dispersion (bottom panel). Dashed lines show the 
relations from our model, equation (Bl), obtained by in- 
serting the values from Tables 1 and 2 in equations (B12) 
and (B9), respectively. The bottom panel can be inter- 
preted as showing that galaxies with large velocity disper- 
sions tend to have older stellar populations. Although the 
steepening at higher redshifts in the bottom panel suggests 
that the scatter around the mean color-magnitude relation 
is larger at high redshift, or that evolution is differential, it 
is, in fact, a selection effect (see text between equations B9 
and BIO). 
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the color-magnitude relation tend to scatter to smaller 
masses from the mass-luminosity relation. The other pan- 
els show that, at fixed mass, redder galaxies tend to be 
fainter (second from top), have smaller sizes (third from 
top), and have larger velocity dispersions (bottom). Since 
luminosity and mass are so tightly correlated, we can think 
of the color-mass residual as being similar to the color- 
magnitude residual. If we treat the color-magnitude resid- 
ual as an age indicator, then Figure 6 indicates that, at 



Fig. 6. — Galaxies which scatter blueward of the color- 
magnitude relation scatter towards smaller masses from 
the mass-luminosty relation (top). At fixed mass, red- 
der galaxies tend to be fainter (second), have smaller 
sizes (third), and have larger velocity dispersions (bot- 
tom). Since mass cx iJ^c^, the trend with a shown in 
the bottom panel could have been inferred directly from 
the trend with size. 
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fixed mass, older (redder) galaxies are fainter, smaller, and 
have larger velocity dispersions. All of these trends are 
qualitatively similar to those of dark matter halos, but 
they are also consistent with a model in which the bluer 
light is less concentrated than the red light (see, e.g. Fig- 
ure 3 in Bernardi et al. 2003a), as color gradients also 
suggest (e.g. Figure 7 in Bernardi et al. 2003d). Figure 7 
is similar, but now the color-magnitude residual is shown 
as the dependent parameter. 

In Figures 6 and 7, the first two panels from the top 
and left deserve further comment. The first panel indi- 
cates that bluer, lower mass galaxies can have the same 
luminosity as redder, more massive galaxies. This is easy 
to understand if the color deviation traces age. Suppose, 
instead, we wanted to construct a model in which the resid- 
uals from the color-magnitude relation trace metallicity 
rather than age, the redder galaxies being the more metal 
rich. In this case, one would interpret the first panel as in- 
dicating that lower mass, metal poor galaxies can have the 
same luminosity as more massive, metal rich galaxies. In 
the Bruzual- Chariot (2003) models, at fixed mass, color 
changes are associated with larger magnitude changes if 
the color change is tied to age rather than metallicity (see 
top- left panel of Figure Al). The two possiblities differ by 
a factor of two, so, in principle, one might to be able to dis- 
tinguish between them using the second panels in Figures 6 
and 7. Figure 6 suggests that {AM\A{g — r)) cx A{g — r), 
and Figure 7 suggests (A(.g - r)\AM) cx AAf/10. These 
observed relations are shallower than both models: shal- 
lower relations indicate that scatter is important, but this 
makes it difficult to determine which of the two possibili- 
ties is more accurate. 

Appendix B develops in some detail the consequences 
of a model in which the assumptions that metallicity in- 
creases with luminosity along the color-magnitude rela- 
tion, and that the residuals from the relation are indi- 
cators of age, are combined with our finding that the 
color-magnitude relation is determined by the color-o" and 
magnitude-cr relations. (The Appendix also shows that 
these conclusions are unchanged if we use the orthogonal 
rather than the direct fit, i.e., if we assume metallicity 
varies along the principal axis of the color-magnitude re- 
lation, and age varies perpendicular to this axis.) In such 
a model: 

• metallicity and a are correlated (galaxies with large 
velocity dispersions tend to be more metal rich); 

• age and a are correlated (galaxies with large velocity 
dispersions tend to be older); 

• lines which run parallel to the mean cr-magnitude 
relation are approximately lines of constant age; 

• in the plane of age along the x-axis and metallicity 
along the j/-axis, lines of constant velocity dispersion 
slope down and to the right. 

If a is indeed sensitive to age, one might wonder if the 
slope of the color-cr correlation steepens with redshift. The 
SDSS is a relatively shallow, magnitude-limited survey; 
the median redshift in our sample is of order 0.1, so the 
sample is not ideal for quantifying trends with redshift. 
On the other hand, it offers a large number of galaxies 



30000) with well-calibrated surface photometry and 
velocity dispersions, so we may reasonably expect to be 
able to detect a small signal. We are also aided by the 
fact that, of the three pairwise relations, color-magnitude, 
(T-magnitude, and color-cr, it is the color-magnitude rela- 
tion which is a consequence of the other two (c.f. Fig- 
ures 2 and 3; also see Bernardi et al. 2003d). Because the 
color-cr relation is the fundamental correlation, selection 
effects associated with the magnitude limit of our sample 
are unimportant (see Appendix B for details). The bot- 
tom left panel of Figure 3 suggests that the slope of the 
color-cr relation is steepening with redshift, but data from 
higher redshift and a better understanding of the appro- 
priate k + e-corrections are required before we can claim 
to have measured this with confidence. Other tantalizing 
evidence of evolution comes from Figure 5: the correlation 
between color-magnitude residuals (i.e., age) and velocity 
dispersion steepens systematically with redshift. Although 
it is tempting to conclude that this is additional evidence 
of evolution, the analysis in Appendix B shows that it is, 
in fact, a selection effect. 

Taken together, the findings itemized above indicate 
that our data require the existence of an age-metallicity 
degeneracy; this degeneracy is most noticable for galaxies 
which have the same velocity dispersion a. In this respect, 
our findings are consistent with Figure 4 of Trager et al. 
(2000b). But our results are inconsistent with their con- 
clusion that the most massive galaxies arc less metal rich 
(see their Figure 6). Indeed, Figure Al indicates that it 
would be difhcult for single burst models to produce the 
observed color-magnitude relation if the more luminous 
galaxies were older but less metal rich. Our finding that 
older galaxies tend to have larger velocity dispersions than 
their luminosities indicate is qualitatively consistent with 
one of the conclusions of Forbes & Ponman (1999); this 
is encouraging because their age estimates come from a 
very different approach (detailed comparison of spectral 
features with stellar population synthesis models). 

While it is tempting to calibrate the relation between 
color-magnitude residual and age, this would be model- 
dependent. For example, the Bruzual-Charlot models pre- 
sented in Appendix A can be used to perform this cali- 
bration. But we are hesitant to take this path because 
the models do not include an important correlation be- 
tween velocity dispersion and a-element abundance ratios 
(Trager et al. 2000a; Thomas, Maraston & Bender 2003; 
Proctor et al. 2004). 

In the model developed here, the scatter around the 
color-magnitude and the a-magnitude relations should in- 
crease with lookback time. This can be tested with data 
from redshifts of order one-half to unity which will soon 
be available (e.g. the SDSS-2dF survey of luminous red 
galaxies at http://sdss2df.ncsa.uiuc.edu/). 
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APPENDIX 



A: SINGLE-BURST MODELS 

In the single- burst models of Bruzual & Chariot (2003), 
galaxy colors depend on both age and metallicity. For 
a population which is older than one Gyr, and has solar 
metallicity or greater, we have found that the color is quite 
well described by: 

g-r = 0.7 + 0.25(r-9.5) + — Z if T > 9.5 

0.40 

g-r = 0.7 + 0.50(r-9.5) + — Z if T < 9.5 



where T — logjQ(age/Gyr). 

The four panels of Figure Al provide different exam- 
ples of the color-magnitude relation associated with these 
models. The top left panel shows the locus traced out 
in the color-magnitude plane by a single object of fixed 
mass as it ages. The two dotted lines which slope from 
top left to bottom right show lines of constant metallic- 
ity (higher metallicity is redder), and the approximately 
vertical lines which connect them show locii of constant 
age (older is fainter). This panel is relevant to results 
presented in Figures 6 and 7. It shows that, at fixed 
mass, color changes are associated with larger magnitude 
changes if the color change is tied to age rather than metal- 
licity: AM (X 5 A{g — r) if the color change is entirely due 
to age, whereas AM oc (5/2) A{g — r) if the color change 
is entirely due to metallicity. 

Top right panel shows models in which metallicity corre- 
lates with luminosity. The solid lines which span approx- 
imately the entire luminosity range show locii of constant 
age in a model which has metallicity increasing with lu- 
minosity. The dashed lines show constant age in a model 
where metallicity decreases with luminosity. Notice that. 
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Fig. A1. — Correlation between color and magnitude in the single-burst models of Bruzual & Chariot (2003). Different 
panels show the result of assuming the various correlations between luminosity, age and metallicity that are discussed 
in the main text. The top right panel shows color magnitude relations in models which have metallicity varying with 
luminosity. The bottom right panel shows relations in which age increases with luminosity. Bottom left panel shows 
models in which age and metallicity vary with luminosity. In the top two panels, solid and dashed lines which connect 
dots trace locii of constant age = 15.25, 12.75, 10.25, 7.75 and 5.25 Gyrs, with the oldest being reddest. 



in both cases, the slope of the color magnitude relation is 
approximately independent of age. However, the data in- 
dicate that color increases with luminosity, so that models 
with metallicity decreasing with luminosity are not per- 
mitted. 

The bottom right panel shows two models in which age, 
rather than metallicity, vary along the luminosity axis. 
Both models assume that older objects are more luminous, 
but the redder set of lines (i.e., the dashed lines) show a 
more metal rich population. In both cases, the slope of 
the relation varies more strongly with age than it did in 
the previous panel, with the slope of the color-magnitude 
relation evolving rapidly when the population is younger, 
but less rapidly at later times. 

The two panels on the right showed models in which 
either metallicity or age varied with luminosity, but not 
both. Dashed lines which slope down and to the right in 
the bottom left panel show the color-magnitude relation 



in a model where the more luminous objects are older but 
metal poor. This model is qualitatively like those of Trager 
et al. (2000b). Solid lines show a model in which both 
age and metallicity increase with luminosity. The data 
are more consistent with this latter model which has more 
luminous objects being older and metal rich. 

B: LUMINOSITIES, COLORS, VELOCITY DISPERSIONS, 
AGES AND METALLICITIES 

Luminosity, color and velocity dispersion 

Define m ~ {M — Ah)/aMM, v = (V — V^,)/avv, and 
c — {C — C^,)/acc, where M is the absolute magnitude, 
V = logjQ (velocity dispersion/km s~^), and C denotes a 
color. Our notation uses X, amd axx to denote the mean 
and rms of observable X = M, V, C, etc. This Appendix 
considers a model in which the only fundamental correla- 
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tions are with a: 

p{v,m,c) — p(v) plmlv) p{c\v) . 



(Bl) 



Pairwise correlations, such as that between luminosity and 
a, are given by integrals like 

(mv) = J dmdvdc{mv)p{v,m,c). (B2) 

For the assumed model, these are relatively simple to eval- 
uate. If we suppose that the mean magnitude at fixed a 
scales as {m\v) = ^mvV, and the mean color at fixed a is 
also linear in v, say ^cvV, then 



{mv) = ^rnv, and {cv) = ^cv 



(B3) 



This comes from our choice to normalize {vv) = 1, and 
similarly for the other variables. 

Let am\v denote the scatter in m around {m\v). If a^iv 
is independent of v, then 



-I 
-J 



dmp(rn) {v'^\m) 



_ 2 , t2 



(B4) 



Therefore, <j'^^^ = 1 — and similarly for the scatter 
around the color-a relation. 

What if we are interested in relations which do not in- 
volve (7? For instance, the mean color at fixed magnitude 
is 



(c|to) 



-I 



dvdccp{c\v) p(v\m) 

dv(,cvV p{v\m) = ^cv^mv m. (B5) 



This shows that the slope of the color-magnitude relation 

can be written in terms of the slopes of the other two 
relations: ^cm = ^.cv^mv This is a simple approximation 
to the observed correlations: ^cm = —0.36 and ^cv^mv = 
-0.40 (from Tables 1 and 2). 

By assumption (Bl), the residuals of the m-a relation, 
Rrn\v = '71— (m|u) = m — S^rnvV, do not Correlate with those 
of the color-cr relation. It is straightforward to verify that 
this is so: 

{Rm\vRc\v) = {{m - ^mvV){c - ^cvV)) 

= {mc - ^cvmv - ^rnvVC + ^mv^cvv'^) 

= 0. (B6) 

However, the residuals of the color-magnitude relations do 
correlate with those of the cr-magnitude relation: 

{Rc\mRv\m) = ((C - ^cTnm){v - (,vm.m)) 

= {cv - ^cmmv - £,yrnmc + ^cmS.vrnm'^) 
~ S,cv ^cmS,mv S,cv (1 S,mv)' i^"^) 

Thus, this model explains the correlations between resid- 
uals shown in Figure 4. 



Assumption (Bl) has another important consequence. 
Namely, the slope of the color-cr relation, at fixed magni- 
tude, is 



f dccp{c\v)p{m\v)p{v) 
{c\v, m) = = ^cvV- 

p[ra, v) 



(B8) 



This shows that the slope of the c-v correlation is the same 
as when the constraint on magnitudes is removed. There- 
fore, the slope of the color-cr correlation measured in a 
magnitude limited sample is the same as the slope of the 
intrinsic relation. In other words, if the only correlation 
between color and magnitude is through cr, then the slope 
of the color-(T relation measured in a magnitude limited 
sample is not biased by selection effects. We use this in 
the main text. 

Age, metallicity and velocity dispersion 

If the residuals from the color- magnitude relation, Rc\rm 
are indicators of age, and if the relation between age and 
Rc\m is linear, then the mean age as a function of velocity 
dispersion can be obtained from 



(i?,|„|t;) = (l-eje 



,v. 



(B9) 



If £,cv > 0, then age increases with increasing velocity dis- 
persion. The rms scatter around this relation is a^v The 

bottom panel of Figure 5 shows this relation in the data. 
In this model, the mean age increases less rapidly with ve- 
locity dispersion (by a factor 1 — ^^^) than does the color, 
suggesting that some of the increase in color is associated 
with factors other than age. 

In Figure 5, the relation between age and a appears 
to steepen with redshift. If real, what does this steep- 
ening imply? Note that the quantities in the Figure 
have not been normalized by their rms. If we suppose 
that the correlations between velocity dispersion, color, 
and magnitude do not evolve, but their rms values do, 
then the steepening implies that {crcc / crvv) is increasing. 
The distribution of velocity dispersions does not appear 
to be evolving (e.g. Sheth et al. 2003), and there are 
no strong reasons to expect it to. On the other hand, 
we argued in the Introduction that one expects to see 
the scatter around the color-magnitude relation increase 
with redshift. Although we do expect to see acc in- 
crease, the increase implied by Figure 5 is larger than 
expected. This is because the steepening is a selection 
effect it is the result of the magnitude limit of the sur- 
vey. To see why, note that the magnitude limit restricts 
the range in m over which the integrals which define 
{Rc\m\''') = Jd,m.p{m,\v) {S,cvV — £,cmfn) arc performed. If 
we write the result as (1 — ^cvV times a correction 
factor / dmp{m\v) {^cvV - ^cmm)/{^ 

CV ^ ^cm ^mvv), then 
this factor depends on the magnitude limit, so in effect, 
the slope of the correlation with v depends on redshift. 

If the residuals from the color-cr relation are a con- 
sequence of age and possibly other variables (such as 
metallicity and/or environment), then correlations be- 
tween Rc\m find Rc\y can be used to study how much of 
the scatter in the color-cr relation is due to age, and how 
much to other effects. For instance, 



{Rc\m\Rc\v} 



R 



'c\v^ 



and 



(BIO) 
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Fig. B1. — Correlations between residuals from the color-magnitude and color-cr relations are well described by the model 
developed in the text. Dashed line in panel on the left shows the correlation described by equation (BIO), and the line in 
the panel on the right shows equation (Bll). 
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(Bll) 



This is consistent with the data has shown in Figure Bl. 
Thus, although objects which have the same fixed offset in 
color from the color-cr relation may have a range of ages, 
the mean age is the same as that infered from the color- 
offset of the color-magnitude diagram (see upper panel of 
Figure Bl). On the other hand, the mean color offset 
from the color-cr relation of objects which have the same 
age (i.e., they are all offset from the color-magnitude rela- 
tion by the same amount) is smaller than the color-offset 
from the color-magnitude relation (sec bottom panel of 
Figure Bl). This is because some the color change asso- 
ciated with age effects has been accounted-for by the fact 
that age varies with a. 

In contrast, the result from Figure 5 (upper panel) show- 
ing 

(i?e|.|m)=0 (B12) 

is less straightforward to interpret. This is because the pre- 
vious correlations suggest that residuals from the color-cr 
relation are some combination of both age and metallicity. 
If luminosity is driven primarily by metallicity, then the 
fact that {R^\y\m) = must mean that the age effects ex- 
actly cancel the correlation with metallicity. Moreover, if 
we replace luminosity with metallicity, then the fact that 
(iJc|ml™) = suggests there is no correlation between age 
and metallicity. 

The mean age at fixed luminosity and velocity dispersion 

is 

{Rc\m\v,m) = ^cv{v - ^mvm), (B13) 

and the rms scatter aromid this mean is ac\v Since the 
final expression is really ^cv Rv\mj this shows that objects 
which have larger velocity dispersions than the mean a- 
magnitudc relation tend to scatter redward of the color- 
magnitude relation, and so are, on average, older than 
those which scatter to smaller dispersions. This has two 
important consequences. First, lines which run parallel to 
the mean {v\m) relation trace out curves of approximately 
constant age, and this statement becomes increasingly ac- 
curate if the relation between color and a is tight (recall 
that the scatter around {Rc\m\v,m) is a^v)- Second, if 
younger objects evolve more rapidly, and the evolution is 
stronger in luminosity than in velocity dispersion, then the 
scatter around the cr-magnitude relation should be larger 
at high redshift. 

The mean velocity dispersion at fixed magnitude and 
color-magnitude residual is 



{v\m, R,\ra) =Uvm+ ^^"/^ .5""^ Rr. 



2 

cm 



(B14) 



with rms scatter cFv\mO'c\v/<^c\m- 

If age and metallicity scale as T = aRc\m and Z = —bm, 
with a and b both positive, then the same average value of 
cr can be obtained by increasing the age and decreasing the 
metallicity, or vice-versa (we have assumed that < 0, 
since, as Figure 1 shows, velocities tend to increase with 
luminosity). We have not discussed the roles played by 
the environment and by changes in a-element abundances. 
In this model, they would be responsible for some of the 
scatter around these mean trends with age and metallicity. 



The expression above suggests that different combina- 
tions of age and metallicity can yield the same average age. 
We can see this more directly by computing the correlation 
between age and metallicity at fixed velocity dispersion: 



{mRc\m\v) - (m\v){R^\^\v) = -C 



(B15) 



If age and metallicity scale as T = ai?c|m and Z = —bm, 
with a and b both positive and £cm < 0, then the same 

value of cr can be obtained by increasing the age and de- 
creasing the metallicity, or vice-versa. This is also consis- 
tent with equation (B13): the mean age at fixed velocity 
dispersion increases as the metallicity decreases. 

So far, we have not asked how one calibrates the relation 
between color-magnitude residual and age. In single burst 
models, it is possible to constrain the ratio a/b, because, 
at fixed cr, login (age) and metallicity are expected to be 
related by a factor of -3/2 (e.g. Worthcy 1994). 

In summary: We showed that the color-magnitude re- 
lation is determined by the color-u and magnitude-cr re- 
lations. Previous work suggests that (i) metallicity in- 
creases with luminosity along the color-magnitude rela- 
tion (as suggested by the observation that the slope of the 
color-magnitude relation does not evolve), and (ii) that the 
residuals from the relation are indicators of age. We de- 
veloped a simple model which combines these results. The 
model implies that: (i) age and a are correlated (galax- 
ies with large velocity dispersions tend to be older); (ii) 
lines which run parallel to the mean cr-magnitude relation 
are approximately lines of constant age; (iii) in the plane 
of age along the x-axis and metallicity along the y-axis, 
lines of constant velocity dispersion slope down and to the 
right; (iv) the scatter around the color-magnitude and the 
(T-magnitude relations should increase with lookback time. 

The orthogonal fit 

The previous section developed a model for correlations 
between age and cr in which luminosity was proportional 
to incitallicity, and residuals from the color-magnitude rela- 
tion were proportional to log(age). Here, we change these 
assumptions slightly, and ask what happens if metallic- 
ity varies along the long axis of the color-magnitude re- 
lation, whereas age varies along the short axis. In this 
case, the age- metallicity plane is simply a rotation of the 
color-magnitude plane, so metallicity is some combination 
of luminosity and color, and age is some different combi- 
nation. In particular, this would set 



Toe 



(JcC - aoo-mm 



and Z oc 



aoacC + amiri 



(B16) 



where 



ao 




When CTc <C <Jm then oq £,cmO'c/o'm, the mean age 
T at fixed velocity dispersion is oc crc(l — Cmt;)Cci;^^5 
and the mean value of Z at fixed u is oc —vam£,mv{'i- + 

^CT^c/'^m)/\/l + ^ -vam£,mv Thus, in this limit (be- 
cause the luminosity in different bands is quite well corre- 
lated, the data do indeed have ac <C cr^), analysis of the 
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orthogonal fit produces results which are similar to the 
analysis of the previous section. 

If the scatter in colors were larger, the difference be- 
tween the orthogonal fit model and the one studied in 
the previous subsection would also be larger. For in- 
stance, if am = cTc, then the age indicator would be 
a (c— m)/-\/2 and the metaUicity would be oc (— c— m)/-\/2, 
so the mean age at fixed velocity dispersion would be 
oc V {S^cv — £,mv)/V^, whereas the mean metallicity at fixed 
V would scale as —v(S,cv + S,mv)/V^- The rms scatter 
around the age-cr and metallicity-cr relations would be 



(cr^i -|- CT^, )/2. The mean correlation between age and 



metallicity at fixed t; is (^^„ — ^^„)/2. In the sample, ^m?; < 
and \^mv\ > S.CV, so the mean age and metallicity increase 
with velocity dispersion, and if the age is larger than ex- 
pected for the velocity dispersion, then the metallicity is 
smaller. Moreover, the mean age would increase as lumi- 
nosity increases: (^cm -l)m/\/2 = -m (1 - £,cv£.mv)/V^, 
with rms scatter ac\rm and the mean metallicity would in- 
crease as {Z\m) = — TO (l-|-^ct)$mt))/'\/2, with the same rms 
scatter, (Tc\m- Although some of these conclusions are qual- 
itatively similar to those of the previous subsection, they 
are generally quantitatively different: the correlations be- 
tween age and velocity dispersion and age and luminosity 
become steeper, whereas those between metallicity and a 
and metallicity and luminosity become shallower. 




